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This  paper  describes  a  solar  heat  pipe  thermoelectric  generator  (SHP-TEG)  unit  comprising  an  evacuated 
double-skin  glass  tube,  a  finned  heat  pipe  and  a  TEG  module.  The  system  takes  the  advantage  of  heat  pipe 
to  convert  the  absorbed  solar  irradiation  to  a  high  heat  flux  to  meet  the  TEG  operating  requirement.  An 
analytical  model  of  the  SHP-TEG  unit  is  presented  for  the  condition  of  constant  solar  irradiation,  which 
may  lead  to  different  performance  characteristics  and  optimal  design  parameters  compared  with  the  con¬ 
dition  of  constant  temperature  difference  usually  dealt  with  in  other  studies.  The  analytical  model  pre¬ 
sents  the  complex  influence  of  basic  parameters  such  as  solar  irradiation,  cooling  water  temperature, 
thermoelement  length  and  cross-section  area  and  number  of  thermoelements,  etc.  on  the  maximum 
power  output  and  conversion  efficiency  of  the  SHP-TEG.  Simulation  based  on  the  analytical  model  has 
been  carried  out  to  study  the  performance  and  design  optimization  of  the  SHP-TEG. 

©  2011  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

With  the  diminishing  reserve  of  fossil  fuels  and  the  pressing  is¬ 
sues  of  environmental  pollution  caused  by  combustion  of  fossil 
fuels,  research  to  find  new  energy  sources  and  affordable  new 
power  generation  methods  has  been  become  more  and  more 
important.  Compared  to  conventional  electrical  power  generator 
systems,  thermoelectric  generators  (TEG)  theoretically  offer  many 
advantages  such  as  being  simple  and  highly  reliable,  having  no 
moving  parts,  and  being  environmentally  friendly.  However,  they 
have  relatively  low  conversion  efficiencies,  so  their  applications 
have  been  usually  limited  to  the  specific  situations  where  reliabil¬ 
ity  is  a  major  consideration  such  as  in  aerospace  and  military  appli¬ 
cations.  More  recently,  there  has  been  a  growing  interest  to  use 
TEG  for  electricity  generation  from  waste  heat  of  various  sources 
such  as  combustion  of  solid  waste,  power  plants,  biomass  cooking, 
and  other  heat-generating  processes  when  the  cost  of  the  thermal 
input  do  not  need  to  be  considered  [1-5],  and  many  researches  on 
system  optimization  and  performance  improvement  of  thermo¬ 
electric  generators  for  waste  heat  recovery  have  been  reported  in 
numerous  publications  [1-10].  In  a  study,  the  cost  of  electricity 
produced  by  TEG  from  waste  warm  water  is  predicted  to  be 
£0.08/kWh  and  even  to  be  £0.04/kWh,  which  can  challenge  the 
cost  of  electricity  produced  by  conventional  methods  from  oil  or 
coal  fuels  [8]. 
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On  the  other  hand,  solar  heat  driven  TEG  is  emerging  as  a  com¬ 
peting  alternative  technology  to  the  dominating  solar  photovoltaic 
(PV)  systems.  Though  its  low  conversion  efficiency  compared  to 
PV,  solar-driven  TEG  technology  still  attracts  increasing  attention 
as  other  simple  and  competitive  way  to  produce  electricity  from 
solar  energy  [11-15].  In  fact,  the  low  conversion  efficiency  may 
be  not  a  serious  barrier  for  use  of  the  free  and  friendly  solar  energy. 
With  the  decreasing  price  of  thermoelement  materials,  the  TEG 
technology  is  attracting  more  research  interest  to  develop  it  for  so¬ 
lar  energy  applications  and  a  number  of  studies  on  design  and  opti¬ 
mization  of  systems  and  analysis  of  their  performance  have  been 
reported  in  the  recent  literature.  The  performance  of  the  TEG  in  a 
hybrid  PV-TEG  system,  where  the  TEG  operates  independently  as 
a  secondary  generator  to  improve  the  overall  efficiency  of  the  sys¬ 
tem,  has  been  studied  [10,11,14].  The  analytical  method  used  in 
these  studies  are  not  much  different  from  that  used  in  the  study 
of  TEG  for  waste  heat  recovery,  and  calculation  of  the  TEG  was 
based  on  the  temperature  difference  between  the  hot  side  and  cool 
side  of  a  TEG,  which  was  usually  chosen  at  60  K,  80  K  or  100  K.  Typ¬ 
ically,  the  thermal  conductivity  of  a  TEG  is  about  1.5  Wm_1  K-1 
and  the  depth  of  commercial  TEG  is  about  4  mm  (including  the 
ceramic  electric  insulation),  so  the  heat  flux  through  a  TEG  could 
be  higher  than  20000  W  m-2  for  a  temperature  difference  of  60- 
100  K.  However,  the  solar  irradiation  is  usually  less  than 
1000  Wm-2,  which  is  obviously  too  low  to  match  the  required 
heat  flux  of  a  TEG  to  obtain  a  large  temperature  difference  for  a 
considerable  efficiency.  Solar  heat  could  be  accumulated  or  solar 
radiation  could  be  concentrated  to  match  the  required  heat  flux 
of  a  TEG.  A  study  has  used  a  liquid  such  as  water  circulation  system 
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Nomenclature 

Symbols 

Subscripts 

A 

area,  m2 

air 

air 

a 

cross-section  area  of  thermoelement,  m2 

c 

convective 

D 

depth  or  diameter,  m 

cop 

plate  copper  fin 

G 

solar  irradiation,  W  m-2 

cross 

the  cross-section  of  plate  copper  fin 

h 

heat  transfer  coefficient,  W  m-2 

load 

external  load 

k 

thermal  conductivity,  W  m_1  K-1 

loss 

heat  loss  from  the  evacuated  glass  tube  to  the  environ¬ 

L 

length,  m 

ment 

l 

length,  m 

in 

internal  load 

P 

power,  W 

r 

radiation 

Q 

heat  flow,  W 

sky 

sky 

R 

load  resistant,  Q 

tube 

glass  tube 

r 

electrical  resistivity,  Q  m 

tube-in 

inner  glass  tube 

T 

temperature 

tube-out 

outer  glass  tube 

U 

voltage,  V 

tegh 

ceramic  plate  of  TEG  hot  side 

W 

width,  m 

tegl 

ceramic  plate  of  TEG  cold  side 

S 

Stefan-Boltzmann  constant 

teg 

TEG 

£ 

emissivity 

f 

water  flow 

T 

transmissivity 

hp 

heat  pipe 

a 

absorptivity  or  Seebeck  coefficient,  V  K_1 

to  transfer  solar  heat  to  a  TEG  while  other  studies  have  investi¬ 
gated  use  of  a  line-focus  parabolic  trough  (PTC),  a  Compound  Par¬ 
abolic  Concentrator  (CPC)  or  a  two-stage  concentrator  to  provide 
high  density  of  solar  radiation  so  as  to  give  a  high  heat  flux  across 
the  TEG  [11,15,16].  A  CPC  usually  has  a  concentration  ratio  of  less 
than  5,  which  is  still  too  low  to  match  the  required  heat  flux  of  a 
TEG,  so  it  needs  to  be  used  along  with  a  circulation  system.  Solar 
tracking  parabolic  concentrators  can  easily  give  a  much  larger  con¬ 
centration  ratio,  but  the  cost  of  a  tracking  system  could  be  a  barrier 
and  also  a  TEG  may  be  not  competitive  to  a  concentrating  PV.  Use 
of  a  liquid  circulation  system  to  transfer  solar  heat  to  a  TEG  is  sim¬ 
ple,  but  a  larger  temperature  change  of  the  liquid  is  required  to  car¬ 
ry  a  certain  amount  of  heat,  and  electricity  consumption  of  the 
circulation  pump  and  heat  loss  of  the  circulation  system  are  other 
concerns. 

And  evacuated  glass  tube  solar  collectors  have  become  common 
practice  in  housing  stock  worldwide,  especially  in  China.  Hot  water 
requirement  in  winter  season  should  be  met  according  to  Chinese 
system  design  principle,  so  it  usually  leads  to  too  much  hot  water 
over  normal  consumption  in  summer  season,  which  leads  to  en¬ 
ergy  waste.  Direct  incorporation  of  a  TEG  with  a  heat  pipe  type 
evacuated  tube  solar  collector  may  offer  another  option  for  design 
of  solar  TEG  systems  and  solar  evacuated  tube  solar  collector  sys¬ 
tems.  The  presented  design  of  a  solar  heat  pipe  TEG  (SHP-TEG)  sys¬ 
tem  is  shown  in  Fig.  1 .  The  SHP-TEG  system  can  supply  hot  water, 
electric  power  or  electric  power  and  hot  water  simultaneously 
according  to  requirement  and  weather  conditions.  This  paper  pre¬ 
sents  design  and  electric  power  performance  analysis  of  a  solar 
heat  pipe  thermoelectric  generator  (SHP-TEG),  which  may  lead  to 
a  high  performance  and  low  cost  solar  TEG  system. 

Compared  with  the  common  study  TEG  which  has  based  the  en¬ 
ergy  balance  and  heat  transfer  analysis  on  any  given  temperatures 
of  the  hot  and  cold  junctions  of  a  TEG,  in  the  study  of  a  SHP-TEG, 
the  solar  irradiation  is  known  as  the  input  to  give  a  constant  heat 
flux,  then  the  hot  and  cold  junction  temperatures  and  other  param¬ 
eters  such  as  the  electric  power  and  efficiency  will  be  determined. 
Hence,  the  thermal  model  and  power  model  in  this  study  will  be 
described  for  the  constant  solar  irradiation  as  one  of  boundary  con¬ 
ditions.  The  effects  of  solar  irradiation,  cooling  water  temperature, 
the  length  and  cross-section  area  of  thermoelements,  and  the  num¬ 
ber  of  thermoelements  on  the  performance  of  a  SHP-TEG  unit  will 


be  analyzed.  This  work  was  intended  to  give  some  guidelines  for 
design  and  operation  of  SHP-TEG  systems. 

2.  Mathematical  models 

It  consists  of  an  evacuated  glass  tube,  a  heat  pipe  with  heat 
transfer  fin  and  a  thermoelectric  module.  The  A-A  cross-sectional 
view  shows  the  details  of  the  thermoelectric  module  and  the  con¬ 
denser  section  of  the  heat  pipe.  A  copper  plate  is  soldered  to  the 
condenser  of  the  heat  pipe  in  order  to  enhance  heat  transfer,  and 
the  thermoelectric  module  may  be  bonded  to  the  copper  plate 
tightly  using  silicon  grease.  The  B-B  cross-sectional  view  shows 
the  detailed  structure  of  the  evacuated  double-skin  glass  tube  with 
the  heat  pipe  being  inserted.  The  evacuated  space  is  between  the 
inner  tube  and  the  outer  glass  tube.  Besides  enhancing  heat  trans¬ 
fer,  the  four  copper  fins  on  the  heat  pipe  are  used  as  the  solar  en¬ 
ergy  absorption  surface  and  also  to  support  the  heat  pipe  inside  the 
inner  glass  tube.  And  the  heat  pipe  converts  low  heat  flux  to  high 
heat  flux  by  changing  the  ratio  of  evaporator  area  to  condenser 
area.  There  are  three  operation  models:  (1)  electrical  power  prior¬ 
ity  model,  the  water  temperature  in  cooling  channel  1  keeps  low, 
the  cooling  channel  2  is  empty  and  closed;  (2)  electrical  power 
and  hot  water  model,  the  water  temperature  in  cooling  channel 
1  meets  requirement,  the  cooling  channel  2  is  empty  and  closed; 
and  (3)  hot  water  model,  thermoelectric  module  does  not  work, 
cooling  channels  1  and  2  are  filled  with  water  and  the  system  oper¬ 
ates  like  a  normal  evacuated  tube-heat  pipe-solar  collector  system. 

The  performance  of  system  operating  under  model  3  like  a  nor¬ 
mal  evacuated  tube-heat  pipe-solar  collector  system  has  been 
studied  in  former  literatures  [17,18].  This  study  is  intended  to  ana¬ 
lyze  the  steady  state  performance  of  the  SHP-TEG  operating  under 
model  1  or  model  2.  The  heat  transfer  process  across  the  SHP-TEG 
is  illustrated  in  Fig.  2.  The  following  assumptions  were  made  to 
simplify  analysis  of  the  complex  heat  transfer  process  in  the 
SHP-TEG: 

(1)  The  fins  of  the  heat  pipe  have  the  same  temperature  as  the 
heat  pipe. 

(2)  All  the  energy  balance  equations  are  for  the  steady  state  con¬ 
ditions,  so  the  thermal  capacity  of  any  parts  of  the  SHP-TEG 
system  is  neglected. 
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A-A 

Fig.  1.  Schematic  of  a  SHP-TEG. 
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Fig.  2.  Schematic  of  the  heat  transfer  network  of  the  SHP-TEG. 


Cooling  water 


(3)  The  contact  thermal  resistance  and  electric  resistance 
between  any  two  parts  of  the  system  are  neglected. 

(4)  The  temperature  gradient  along  the  heat  pipe  is  neglected. 

(5)  The  section  containing  the  condenser  of  the  heat  pipe,  the 
TEG  and  the  water  jacket  are  well  insulated  without  any 
heat  loss  to  the  surroundings. 

2.1.  Evacuated  glass  tube 

The  heat  balance  of  the  evacuated  double-skin  glass  tube  may 
be  given  by: 

Q-Zoss  ~~  hC,air(TtUbe  —  T air)  +  hrxky(Ttube  —  Tsky) 

—  hC,tUbe-hp(Thp  —  Tfube)  +  hr,tUbe-hp(Thp  —  T  tube)  (1) 

where  QiOSs  is  the  heat  loss  from  the  evacuated  glass  tube  to  the 
environment,  Ttube  is  the  surface  temperature  of  the  outer  tube, 
Thp  is  the  temperature  of  the  heat  pipe,  fins  and  the  inner  glass  tube, 
Tair  is  the  temperature  of  ambient  air,  Tsky  is  the  sky  temperature, 
hca  is  the  convective  heat  transfer  coefficient  between  the  surface 
of  the  outer  glass  tube  and  the  ambient,  hrxky  is  the  radiative  heat 
transfer  coefficient  between  the  surface  of  the  outer  glass  tube 
and  the  sky,  hcXube.bp  is  the  convective  heat  transfer  coefficient  be¬ 
tween  the  outer  and  inner  glass  tubes,  which  can  be  assumed  to 


be  zero  because  of  the  evacuated  space,  and  hrXube.bv  is  the  radiative 
heat  transfer  coefficient  between  the  outer  and  inner  glass  tube. 

The  convective  heat  transfer  coefficient  between  the  surface  of 
the  outer  glass  tube  and  the  ambient  may  be  given  by 
hc, air  =  3.8  +  2. 7  wind. 

The  radiative  heat  transfer  coefficient  between  the  surface  of 
the  outer  glass  tube  and  the  sky  may  be  given  by  hrxky  =  stubeS 

(jsky  +  (j sky  +  Tfube')- 

The  convective  heat  transfer  coefficient  between  the  outer  and 
inner  glass  tubes  may  be  given  by: 

hr,tube-bp  —  ^(jtube  +  i^tube 

+  Thp)/(  1  /  Wm  +  (1  /  Wout  -  1 )  j 

where  S  is  the  Stefan-Boltzmann  constant  and  stube  and  sbp  are  the 
emissivity  of  the  outer  glass  tube  and  inner  glass  tube,  respectively, 
Atube-out  is  the  surface  area  of  the  outer  glass  tube,  Atube_in  is  the  sur¬ 
face  area  of  the  inner  glass  tube.  There  are: 

Atube-out  —  ^Dtube-out^tube  cind  Afube-in  =  tube-in\-tube 

where  Dtube.out,  Dtube.in ,  Ltube  are  the  diameters  of  the  outer  and  inner 
tubes  and  the  length  of  the  evacuated  glass  tube  respectively. 
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2.2.  Heat  pipe 

The  temperature  gradient  along  the  heat  pipe  and  its  fins  is  ne¬ 
glected.  The  solar  energy  absorbed  by  the  heat  pipe  is  converted 
into  heat,  most  of  which  is  transferred  to  the  condenser  part  of 
the  heat  pipe  where  the  TEG  high  temperature  side  is  attached  to 
and  the  rest  is  transferred  to  the  outer  glass  tube  and  then  lost 
to  the  ambient.  The  energy  balance  of  the  heat  pipe  may  be  written 
as: 


Qteg-f  —  htegl-f^tegl  ( T  tegl  Tf)  —  2  ElfegCCtegT  tegJ 

+  2nteg  C- ^  ( Ttegh  -  Tteg:)  +  h22 nteg  r~^l  (7) 

l  teg  ^  tit  eg 

where  htegi-f  is  the  effective  heat  transfer  coefficient  between  the 
ceramic  plate  and  the  cooling  water  and  it  may  be  set  to  be  an 
equivalent  constant  value,  for  example  2000  W  m-2  K-1  due  to 
water  cooling  and  heat  transfer  enhancement  of  fins,  and  7)  is  the 
temperature  of  cooling  water. 


GT  mbe-out  Ottube-inAhp- a 


hr,tube-hpAhp-oc  (' Thp  ~  T  tube ) 
"T  h-hp-teghAcoph  (Thp  ~  Ttegh) 


(2) 


where  G  is  the  solar  irradiation  density,  t tube-out  is  the  transitivity  of 
the  outer  glass  tube,  oitUbe-in  is  the  absorptivity  of  the  inner  glass 
tube,  Ahp_a  is  the  heat  pipe  absorber  area,  Ahp_a  =  Dtube.inLtube ,  Acoph 
is  the  plate  copper  fin  area,  Acop  =  WcopLcop,  Wcop,  Lcop  is  the  width 
and  length  of  plate  copper  fin,  Ttegh  is  the  temperature  of  ceramic 
plate  (assumed  to  be  same  as  the  temperature  of  TEG  hot  side), 
and  the  width  and  length  of  ceramic  plate  are  same  as  that  of  plate 
copper  fin.  hhp.tegh  is  the  heat  transfer  coefficient  between  the  heat 
pipe  condenser  and  the  plate  ceramic  plate  of  the  TEG,  as  shown  in 
Fig.  1B-B.  It  is  a  tube-fin  structure  and  the  end  of  the  plate  copper 
fin  is  well  insulated,  so  the  heat  transfer  coefficient  can  be  written 
as  following: 


h-hp-tegh  —  2 y  hCop-teghTTcopkCopAcross  tanh 

/WPlcoplcop)  /  {Acop/2)  (3) 

Itcopficross 

where  hcop-tegh  is  the  heat  transfer  coefficient  between  the  plate  cop¬ 
per  fin  and  the  ceramic  plate  and  it  is  just  the  average  thermal  con¬ 
ductivity  as  expressed  as  following: 

hcop-tegh  —  1/(2  DC0p/kC0p  +  2D  tegh  /  ktegh)  (4) 

where  Dcop  and  kcop  is  the  depth  and  conductivity  of  the  copper  fin 
respectively,  Dtegh  and  ktegh  is  the  depth  and  conductivity  of  the 
ceramic  plate  respectively,  and  Across  is  the  cross-section  area  of 
plate  copper  fin,  Across  =  VVC0pT)C0p,  PLcop  =  2 (WCOp  +  Dcop )  is  the  perim¬ 
eter  of  the  cross-section  of  the  plate  copper  fin  and  it  may  be  sim¬ 
plified  to  be  PLcop  =  2Wcop  neglecting  Dcopt  then  Eq.  (3)  can  be 
rewritten  as: 


hhp-tegh  —  2 y  2hC0p-teghkC0pDC0p  tanh 

X  (^W( 2k^D^)!^j /(Lav/ 2)  (5) 


2.3.  The  hot  side  of  TEG 

Heat  balance  of  the  hot  side  of  the  TEG  may  be  given  by: 

< hr  -teg  —  hhp-tegh  Acoph  {T hp  Ttegh)  —  2  El  teg  OC  teg  T  tegh  I 

+  2nteg [Ttegh  ~  Ttegl)  - \l22 nteg  (6) 

l  teg  2  a  teg 

where  nteg  is  the  number  of  thermoelements,  ateg  is  the  Seebeck 
coefficient  of  thermoelectric  materials  (p-type  and  n-type  assumed 
to  be  same),  /  is  the  current,  ateg  and  lteg  are  the  cross-section  area 
and  length  of  thermoelement  respectively,  rteg  and  kteg  are  the  elec¬ 
trical  resistivity  and  thermal  conductivity  of  thermoelectric  materi¬ 
als,  and  Ttegi  is  the  temperature  of  the  TEG  cold  side). 

2.4.  The  cold  side  of  TEG 


2.5.  Electrical  voltage  and  power 

The  electrical  open  circuit  voltage  of  the  TEG  is  given  by: 

D  =  2 ritegOCteg (Ttegh  ~  Ttegl)  (8) 

The  electrical  power  output  is: 

P  =  U  Pload/ (Pload  +  2 EltegE  tegheg  /  Q-teg)  (9) 

Rioad  is  the  load  resistance.  The  electrical  power  efficiency  of  the 
SHP-TEG  may  be  defined  as  the  ratio  of  the  electrical  power  output 
to  the  received  solar  radiation,  that  is: 

t]  —  P  /  ( GAtube  -out  )  (10) 

3.  Simulation 

The  mathematical  model  described  above  was  written  into  a 
FORTRAN  programme  to  simulate  the  performance  of  the  SHP- 
TEG  system  for  different  design  and  operating  conditions.  The  geo¬ 
metrical  and  physical  parameters  of  the  TEG,  the  heat  pipe  and  the 
evacuated  double-skin  glass  tube  were  chosen  referring  to  the  gen¬ 
eral  specifications  of  commercially  available  components.  The 
parameters  used  in  simulation  are  listed  in  Tables  1-3  referring 
to  commercial  production  and  those  physical  parameters  were  as¬ 
sumed  to  be  constant  in  the  chosen  range  of  temperatures.  The  TEG 
was  comprised  of  a  number  of  Bi2Te3  thermoelements  with  differ¬ 
ent  cross-section  areas.  The  basic  conditions  of  simulation  were  set 
as:  the  solar  irradiation  is  1000  W  m-2,  the  cross-section  area  and 
the  number  of  thermoelements  are  1.4  x  1.4  mm  and  127  respec¬ 
tively,  the  ambient  temperature  is  25  °C,  the  wind  velocity  is 
3  ms"1,  and  the  TEG  cooling  fluid  is  water. 

For  the  condition  of  a  constant  AT,  the  previous  studies  in  the 
literature  have  reported  that  the  maximum  electrical  power  of 
TEG  can  be  achieved  when  the  load  resistance  equals  the  internal 


Table  1 

Geometrical  and  physical  of  the  TE  model. 


CCteg 

l<teg 

Vteg  111)  ktegh 

El  tegh 

El  tegl 

(VK-1) 

(W  ITT1  K1) 

(W  ITT1  K1) 

(m) 

(m) 

2.0  x  10“4 

1.5 

9.5  X  10“6  25 

0.001 

0.001 

Table  2 

Geometrical  and  physical  parameters  of  the  heat  pipe. 


thp  (m) 

Wcoph  (m) 

Ecoph  (m) 

Dcoph  (m) 

kcop  (W  m~’  K~’) 

1.8 

0.04 

0.04 

0.002 

398 

Table  3 

Geometrical  and  physical  parameters  of  the  evacuated  glass  tube. 

Etube  (m)  Dtube  -out  (m)  El^ube  -in  (m)  £ tube-out  ^ tube-out  & tube-in  & tube-in 

1.8  0.057  0.047  0.88  0.95  0.06  0.95 


Heat  balance  of  the  cold  side  of  TEG  may  be  given  by: 
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Internal  Resistance  (ohm) 

Fig.  3.  External  load  matches  internal  resistance  to  get  the  maximum  power  output 
for  various  solar  irradiation  values  ( nteg  =  127,  7)=  298  K). 


Fig.  4.  AT  vary  with  load  resistance  for  various  solar  irradiation  values 
(lteg  =  2.5  mm,  nteg  =  127,  Tf=  298  K). 


electrical  resistance.  But  for  the  condition  of  constant  solar  irradi¬ 
ation  dealt  with  in  the  current  study,  the  load  resistance  Rioad  is  not 
the  same  as  the  internal  electrical  resistance  rin  when  the  maxi¬ 
mum  electrical  power  is  achieved,  but  a  little  larger  than  the  inter¬ 
nal  electrical  resistance,  as  shown  in  Fig.  3.  This  could  be 
demonstrated  through  the  following  deduction. 

Substituting  Eq.  (8)  into  Eq.  (9)  gives  a  new  expression  of  the 
power  output  of  a  TEG: 

n  Up  (2  ritegOCteg)  AT 

( Rioad  +  Rin )  ( Rioad  +  Rin ) 

where  the  internal  electrical  resistance  is: 

Rin  —  2  H-tegT  tegheg  /  Q-teg 

Differentiating  Eq.  (11)  gives: 

dP  (2 TltegOC^eg)  Rioad  d(AT  )  (2 HtegO^teg)  AT 

dRload  (Rioad  +  Rin)2  dRioad  (Rioad  +  Rin)2 
2  (2 Tlteg^teg)  AT  Rioad 

(Rioad  +  Rin)3 


(11) 

(12) 

(13) 


which  should  be  equal  to  zero  in  order  to  achieve  the  maximum 
electrical  power  output. 

For  a  constant  AT,  if 


dP  =  (2nteguteg)2 AT2  Ri"  R,Mi  T  =  0 


dRload 

then, 


(Rioad  +  Rin) 


Rioad  —  Rin 


(14) 


(15) 


For  the  condition  of  constant  solar  irradiation  G,  AT  varies  with 
Rioad  and  the  slope  is  positive,  e.g.  >  0,  as  shown  in  Fig.  4,  so  if 

dP  2  AT2 Rin  —  AT2 Rioad  +  2 Rioad  (Rioad  +  Rin)  AT 

Tfp - “  (2nteg0Cteg)  - — - ;  -  ~  0 

dRload  (Rioad  +  Rin) 

(16) 


Then, 

2dAT  2 

A TdRioad  load 


T2RindAT 

\ATdRioad 


Rioad  +  Rin  —  0 


(17) 


which  is  complex  to  get  Rioad  expression  with  Rin  as  variable.  In  or¬ 
der  to  get  a  simple  expression  of  R[oad ,  R2oad  was  assumed  to  be  RinR- 
ioad •  As  there  is  not  big  difference  between  Rin  and  Rioadt  this 


Fig.  5.  Max  power  efficiency  changes  with  TEG  length  for  various  solar  irradiation 
values  (ambient  temperature  =  25  °C  and  cooling  water  temperature  =  25  °C). 


assumption  would  not  cause  significant  error,  so  Eq.  (17)  can  be 
written  as: 


2  dAT  D  D  ,  (2 Rin  dAT 
ATdRload  m  load  +\AT  dRload 


Rioad  A  Rin  —  0 


(18) 


Thus.  R,00d  =  i  _^j^>Rm  (19) 

AT  dRioad 

Eq.  (19)  shows  a  clear  relationship  among  Riw  Rioad ,  AT  and  As 

UKload 

shown  in  Fig.  4,  the  does  not  change  much  with  variation  of  the 
solar  irradiation.  According  to  Eq.  (19),  when  a  lower  solar  irradia¬ 
tion  leads  to  a  lower  AT,  a  larger  matched  load  resistance  will  be 
needed,  as  seen  in  Fig.  3. 

Fig.  5  shows  variation  of  the  SHP-TEG  power  conversion  effi¬ 
ciency  with  the  change  of  the  solar  irradiation  and  thermoelement 
length.  It  can  be  seen  in  Fig.  5  that  the  maximum  power  conversion 
efficiency  increases  with  increasing  solar  irradiation  density.  For 
the  condition  of  constant  temperature  difference  between  the 
hot  and  cold  sides  of  a  TEG,  the  conversion  efficiency  usually  in¬ 
creases  with  increasing  the  thermoelement  length  and  the  maxi¬ 
mum  power  output  will  be  reduced.  But  for  the  condition  of 
constant  solar  irradiation  discussed  here  the  maximum  conversion 
power  efficiency  has  a  peak  point  when  varying  the  thermoele¬ 
ment  length  and  so  does  the  maximum  power  output  according 
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Fig.  6.  AT  and  Qi0SS  vary  with  the  thermoelement  length  for  various  solar  irradiation 
values  and  rin  vary  with  the  thermoelement  length. 

to  Eq.  (10).  For  a  given  solar  irradiation  which  causes  a  constant 
heat  flux,  a  longer  thermoelement  length  means  a  larger  thermal 
resistance  and  a  larger  AT,  and  thus  leads  a  larger  operation  volt¬ 
age  and  higher  conversion  efficiency.  However,  when  AT  is  larger, 
the  temperature  of  hot  junction  also  becomes  higher,  so  the  oper¬ 
ation  temperature  of  heat  pipe  is  higher  and  heat  loss  Q/oss  from  the 
heat  pipe  to  the  ambient  becomes  larger.  When  the  length  of  ther¬ 
moelement  reaches  a  certain  value,  the  effect  of  heat  loss  will  be¬ 
come  larger  than  the  effect  of  increased  AT,  and  then  the  power 
efficiency  becomes  to  decrease  with  the  length  of  thermoelement. 
This  could  be  further  explained  referring  to  Fig.  6  which  shows  AT 
and  Q/oss  variation  with  the  length  of  thermoelement  for  different 
solar  irradiation.  For  example,  under  the  solar  irradiation  of 
1000  W  m-2,  the  slope  of  the  curve  of  Q/oss  is  lower  than  that  of 
the  temperature  difference  when  the  thermoelement  length  is  less 
than  2.5  mm.  The  two  slopes  are  almost  the  same  when  the  ther¬ 
moelement  length  is  around  2.5  mm  and  then  begin  to  separate 
after  2.5  mm.  This  indicates  that  the  thermoelement  length 
2.5  mm  could  be  a  transition  point  to  the  performance  characteris¬ 
tics  of  the  SHP-TEG  and  might  be  the  reason  why  the  peak  point  of 
conversion  efficiency  appears  at  the  thermoelement  length  of 
about  2.5  mm.  under  solar  irradiation  1000  W  m-2.  It  is  also  evi¬ 
dent  from  Fig.  5  that  the  position  of  the  peak  point  of  conversion 
efficiency  will  change  if  the  solar  irradiation  is  different.  For  exam¬ 
ple,  under  the  solar  irradiation  of  600  W  m-2,  the  peak  point  ap¬ 
pears  at  the  thermoelement  length  of  about  3.5  mm  compared 
with  2.5  mm  for  1000  W  m-2.  A  smaller  solar  irradiation  will  give 
a  lower  heat  flux  across  the  TEG  and  thus  a  smaller  AT  and  lower 
operation  voltage.  As  a  larger  thermoelement  length  can  lead  to  a 
larger  operation  voltage,  the  thermoelement  length  will  need  to  be 
longer  to  achieve  large  enough  operation  voltage  under  the  lower 
solar  radiation. 

The  internal  electrical  resistance  and  thermal  resistance  of  TEG 
also  depend  on  the  number  of  thermoelements  and  their  cross-sec¬ 
tion  area,  and  the  electric  power  output  will  vary  with  the  internal 
electric  resistance  and  thermal  resistance  of  TEG  as  analyzed  above. 
Fig.  7  shows  the  maximum  power  efficiency  of  the  SHP-TEG  varies 
with  the  number  of  thermoelements  for  various  cross-section  areas 
of  thermoelements  under  the  conditions  of  solar  irradiation 
1000  W  m-2  and  thermoelement  length  25  mm.  As  seen  in  Fig.  7 
that  the  optimal  number  of  thermoelements  decreases  with  the 
cross-section  area  increasing.  However,  regardless  of  how  the  num¬ 
ber  of  thermoelements  and  cross-section  area  change,  the  maxi¬ 
mum  power  efficiency  they  can  achieve  is  the  same  and  it  is 
3.346%  for  the  simulation  conditions.  The  optimal  number  of 


Fig.  7.  Maximum  power  efficiency  varies  with  the  number  of  thermoelements  for 
various  cross-section  areas  of  thermoelement  (G  =  1000  W  m-2,  !teg  =  2.5  mm, 
Tf=  298  K). 


thermoelements  is  139  for  1.3  mm  x  1.3  mm,  120  for  1.4  x  1.4, 
105  for  1.5  x  1.5,  and  92  for  1.6  x  1.6,  Multiplying  the  optimal 
number  of  thermoelements  and  their  cross-section  area  gives  the 
same  total  area  of  thermoelements,  i.e.,  about  235  in  this  case.  This 
means  that  to  achieve  the  same  maximum  electrical  power  effi¬ 
ciency  the  same  amount  of  thermoelement  materials  would  be  re¬ 
quired.  This  provides  an  important  design  information  for  the  cost 
calculation  of  the  SHP-TEG  system  as  the  cost  of  thermoelement 
materials  is  a  major  part  in  the  total  cost  of  a  SHP-TEG  system. 

Fig.  8  shows  the  maximum  power  efficiency,  the  hot  and  cold 
side  temperatures  of  TEG  and  the  temperature  of  heat  pipe  vary 
with  the  temperature  of  cooling  water  for  G=  1000Wm-2, 
lteg=  2.5  mm,  nteg  =  120.  The  maximum  power  efficiency  decreases 
from  3.346%  to  2.611%  and  the  temperature  of  the  cold  side  of  TEG 
increases  about  30  K  when  the  temperature  of  cooling  water  in¬ 
creases  from  to  298  K  to  328  K,  and  meanwhile  the  hot  side  tem¬ 
perature  of  TEG  is  increased  by  about  10I<  which  is  not  as  large 
as  that  of  the  cold  side  of  TEG.  The  maximum  power  efficiency  de¬ 
creases  about  0.12%  for  every  5  K  increase  in  the  temperature  of 
cooling  water,  so  the  balance  needs  to  be  considered  carefully 
when  the  SHP-TEG  is  designed  for  combined  electricity  and  heat 
supply  applications.  Furthermore,  when  the  temperature  of  cooling 
water  is  328  K,  the  temperature  of  the  hot  side  of  TEG  reaches 
450  I<  which  nearly  touches  the  melting  point  456  K  of  tin/lead  sol¬ 
der  and  risks  damaging  the  TEG. 


Fig.  8.  Maximum  power  efficiency  and  temperatures  of  TEG  vary  with  temperature 
of  cooling  water  (G  =  1000  W  m-2,  lteg  =  2.5  mm,  nteg  =  120). 
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4.  Conclusions 

An  analytical  model  of  the  solar  heat  pipe  thermoelectric  gener¬ 
ator  (SHP-TEG)  system  has  been  developed  for  the  condition  of 
constant  solar  irradiation  based  on  operation  models  1  and  2, 
which  causes  difference  performance  characteristics  and  optimal 
design  parameters  compared  with  the  condition  of  a  constant  tem¬ 
perature  difference  across  a  TEG.  The  analytical  model  presents  the 
complex  influence  of  basic  parameters  such  as  solar  irradiation, 
cooling  water  temperature,  thermoelement  length  and  cross-sec¬ 
tion  area  and  number  of  thermoelements,  etc.  on  the  maximum 
power  output  and  conversion  efficiency.  The  analytical  model  of 
the  SHP-TEG  provides  a  tool  for  design  optimization  and  perfor¬ 
mance  prediction  for  the  SHP-TEG  system.  Based  on  the  analytical 
model,  simulation  has  been  carried  out  to  study  the  performance 
and  design  of  the  SHP-TEG.  The  simulation  shows  that  there  exists 
a  peak  point  of  the  maximum  power  conversion  efficiency  when 
the  thermoelement  length  changes  and  the  efficiency  of  3.346% 
can  be  achieved  for  1000  W  m-2  solar  irradiation.  To  get  the  max¬ 
imum  power  efficiency  in  the  condition  of  constant  solar  irradia¬ 
tion,  the  matched  external  load  resistance  is  a  little  larger  than 
the  TEG  internal  electrical  resistance.  This  is  difference  from  the 
condition  of  constant  temperature  difference  across  a  TEG,  where 
the  matched  external  load  resistance  should  be  equal  to  the  inter¬ 
nal  electrical  resistance  of  a  TEG. 

It  is  worth  mentioning  that  the  factors  associated  with  manu¬ 
facturing  and  fabrication  processes  such  as  contact  thermal  resis¬ 
tance,  contact  electric  resistance  and  temperature  dependence  of 
material  properties,  degrade  of  material  properties,  etc.  have  been 
neglected  in  this  study,  but  they  are  also  important  to  the  SHP-TEG 
system  and  could  be  considered  in  the  further  performance  study 
of  the  SHP-TEG. 

Furthermore,  the  transient  state  performance,  long  term  perfor¬ 
mance  considering  degrade  of  material  properties  and  annual  per¬ 
formance  analysis  of  the  SHP-TEG  system  with  switching  operation 
model  will  be  performed  in  next  step.  As  there  is  a  big  market  of 
evacuated  glass  tube  solar  collector  in  China,  and  it  is  easy  to  man¬ 
ufacture  the  SHP-TEG  system  technically  based  on  a  solar  evacu¬ 
ated  tube  solar  collector  system,  with  the  decreasing  price  of 
thermoelement  materials,  SHP-TEG  system  will  be  economical  in 
the  near  future. 
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